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Abstract. Dipole, triangular, and higher harmonic flow that have an origin in the 
initial density fluctuations has gained a lot of attention as they can provide additional 
important information about the dynamical properties (e.g. viscosity) of the system. 
The fluctuations in the initial geometry should be also reflected in the detail shape and 
velocity field of the system at freeze-out. In this talk I discuss the possibility to measure 
such fluctuations by means of identical and non-identical particle interferometry. 

The initial spatial conditions in heavy ion collisions are not azimuthally symmetric. 
Due to particle rescatterings it leads to the anisotropic flow - anisotropics in particle 
momentum distributions - with elliptic flow being the well known example. Recently, 
a significant progress has been reached in understanding the role of the initial density 
fluctuations. In particular it was realized that such fluctuations lead to odd harmonic 
anisotropic flow, studying of which brings new insights into dynamics of the system 
evolution. Below, I briefly summarize recent developments. I apologize for not providing 
the full list of references - many of which can be found in other proceedings of this 
conference. 

An importance of the initial state fluctuations for flow development has been noticed 
in PQ, though the exact relation between fluctuation and the final event anisotropics 
were not clear. In [2] an observation was made that non-zero two-particle rapidity 
correlations (as observed in pp collisions) in conjunction with radial flow (accepted part 
of the system evolution in AA collisions) lead to a narrow in azimuth and long ranged in 
rapidity correlations. Such correlations were observed experimentally and called ridge. 
This mechanism (non-zero "primordial" rapidity correlations + radial flow =4> ridge) 
have been exploited later in several other models that try to address in more detail a 
(more subtle) question of the origin of the "primordial" rapidity correlations. Although 
the appearance of the ridge and anisotropic flow fluctuations look as totally unrelated 
phenomena, they appeared to be different views on the same thing - reaction of the 
system to the fluctuation in the initial density. In [3j it was shown explicitly that the 
fluctuations in the initial density distribution that extends over large rapidity range 
lead to the ridge structure in two particle correlations. It was also noticed that after 
subtraction of the expected contribution from elliptic flow, two-bump structure appears 
on the away-side. Further studies [4J with a single "hot spot" embedded into otherwise 
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smooth background appeared quite revealing. Contrary to an expectation [2] that such 
a hot spot moved out by the radial flow would create a "bump" in azimuthal distribution 
it appeared that that the hot spot actually "blocks" the development of the radial flow 
and leads to a "dip" in particle distribution at the corresponding azimuth, accompanied 
by side-splashes from both sides. Remarkably, in terms of the correlation function the 
"bump" and the "dip" leads to the same structure - the ridge, as both means positive 
correlations. At the same time, the exact shape of the correlation function (e.g. the 
strength of the third harmonic, which is important below) is different. 

It was noticed in [6] that the fluctuating initial conditions generate anisotropic 
flow of different harmonics. That followed by understanding El E] that the 
perturbations due to each of the hot spots can be treated independently, which allows to 
reformulate the problem from a different perspective - decompose the initial density into 
multipoles and study the system responce to each of the multipoles - the approach widely 
discussed at this conference as flow (dipole, triangular, quadrangular, etc.) fluctuations. 
One can envision such an approach as rotating each of the event to a given harmonic 
symmetry plane, such that at the end one could have a smooth initial conditions but 
with a shape corresponding to a given harmonic. The question addressed in this talk is 
if with the help of femtoscopy (identical and non-identical two-particle correlations) one 
can directly observe such triangular etc. shapes of the system. In this study I follow an 
approach of the first papers where the azimuthally sensitive femtoscopic analyses have 
been proposed [HUE], and blast wave and AMPT models calculations. 
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Figure 1. Left: side-out coordinates. Middle: illustration of predominant expansion 
along shorter directions. v n (pt) for typical values of parameters used in this work. 

Intensity interferometry (identical particle femtoscopy), measures particle space- 
time distribution at freeze-out. Due to the probe on-mass constraint, the full space-time 
information can not be obtained, and one can access only the distribution in (r — Vi), 
where r is the particle position at time t (after freeze-out), and V is the particle velocity. 
For a Gaussian source, the correlation function appears to be a Gaussian 

C(q,P)cxl + expXX-fcfc, (1) 

h3 

where q = pi — p 2 is the pair relative momentum, P = (pi + p 2 )/2 ~ p^ is the 
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particle average momentum, and i??- = ((r^ — Vit)(rj — Vjt)) are the HBT radii. The 
deviation from a Gaussian form can be studied by evaluating the higher moments of 
the distribution (for the effect of non-Gaussiness on HBT radii, see [12] ). Detailes of 
femtoscopic analyses can be found in a review [13J. 



harm 2, a 


e =0.05, a boos =0.05 


— < x L> " 


p = [0.2,0.4] GeV 


— <*L>-<x o 


t > 2 .... p = [0.4,0.6] GeV 




p = [0.6,0.8] GeV 




0.5 



1.5 



2.5 



3 



£ 50 

CM ■=• 

IQI 40 
30 
20 
10 


-10, 





1 1 1 1 1 1 1 


, | , , , | , , , | , , , | , , , | , , , | , , , | , , , 






harm 3 


, a, =0.02, a =0.03 

space boost 






(X 2 

V Sid 


e } p = [0.2,0.4] GeV 






— (x 2 u 


>.(x out > 2 ---- p = [0.4,0.6] GeV 






— (\« 


. * out r P = [0-6,0.8] GeV 






























~l 1 1 I I 1 I I I 1 I I 1 1 1 I 1 I 1 I 1 I 1 I 1 I 1 I I 1 I I I 1 I I I 1 1 I l~ 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



£ 50 

cm 

WH 40 

30 
20 
10 


-10, 



harm 4, a space =0.02, 


a boos =0.03 




p = [0.2,0.4] GeV 




p = [0.4,0.6] GeV 


— < X sid e X -t) 


p= [0.6,0.8] GeV 



0.2 0.4 0.6 0.8 1 1.2 1.4 



Figure 2. Blast wave model: HBT radii azimuthal dependence relative to the 
corresponding symmetry planes for harmonics 2, 3, and 4. 

We use a standard side-out-long system [13J (see figure [j]) and concentrate on 
azimuthal dependence of the side, out, and side-out radii. For a stationary (not 
expanding) source the radii azimuthal dependence can be expressed as 

Ride = inside) = (x 2 ) sin 2 <f> + (y 2 ) cos 2 <j> - (xy) sin 20, (2) 
which has only n = 2 harmonic. Higher harmonics azimuthal dependence appears 
only as deviation from the Gaussian, e.g. in (x® ide ) and (x% de ) for triangular and 
quadrangular shapes respectively. The picture changes if one considers azimuthal 
variation in the expansion velocity. To study this effect we employ a blast wave 
model. We assume longitudinally boost invariant source and freeze-out at a constant 
temperature T = 100 MeV. Expansion velocity profile is parametrized in a form 

T 

pt(r, (f)) = p t ,max 7tt[1 + a boost cos (n(f))} , (3) 

rmax = R[l + a space cos(n0)] , (4) 
with direction of the collective velocity being perpendicular to the line oc [1 + 
a S pace cos(n0)] going through the emission point. The parameter pt, max = 0.9 which 
corresponds to an average expansion velocity (v r ) « 0.7 and (pt)^ = 0.40 GeV. a space 
and ciboost parameters define the spatial eccentricity at freeze-out and modulation in 
the expansion velocity. For the first we assume that the system spatial eccentricity is 
reduced about two times during the system evolution, which leads to a space ~ 0.05; ab oost 
can be estimated from v n (p t ) dependence and is taken to be in the range 0.03-0.07. The 
azimuthal dependence of the obtained HBT radii, see figure [2j indicates that the higher 
harmonics shape effects become clearly visible (and measurable). Note the difference 
in the amplitudes of R s ide-out terms. The radii modulation strongly depends on the 
parameters of the model, and the correspondng measurements will be very important 
in evaluating the velocity profiles at freeze-out and testing models. 
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The AMPT model [14] has been extensively used in studying the effect of the initial 
state fluctuations. We use this model to study Au+Au collisions at ^/5 NN =200 GeV. 
Figure [3] shows the HBT radii variation as function of the azimuth relative to a given 
harmonic event plane which includes the effects of finite reaction plane resolution. Once 
again, the radii dependence on the azimuth is clearly seen. A detailed investigation 
of the side-out relative phases, dependece on the particle transverse momentum, etc. 
requires additional investigations. 
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Figure 3. AMPT model: the azimuthal dependence of the HBT radii (left and 
middle), and radial shift in the production points of kaons and pions, x out = x out — V t t, 
at p t ~ 1.5 rn. 



Summary. Using the blast wave and AMPT models we have demostrated the possibility 
of femtoscopic higher harmonics system shape analysis that promices new insights into 
the physics of heavy ion collisions. 

References 

[1] C. E. Aguiar, Y. Hama, T. Kodama, T. Osada, Nucl. Phys. A698, 639-642 (2002). 

[2] S. A. Voloshin, Phys. Lett. B632, 490-494 (2006). Nucl. Phys. A749, 287-290 (2005). 

[3] J. Takahashi et al. Phys. Rev. Lett. 103, 242301 (2009). 

[4] R. P. G. Andrade, F. Grassi, Y. Hama, W. -L. Qian, Nucl. Phys. A854, 81-88 (2011). 

[5] A. P. Mishra, R. K. Mohapatra, P. S. Saumia, A. M. Srivastava, Phys. Rev. C81, 034903 (2010). 

[6] P. Sorensen, J. Phys. G G37, 094011 (2010). 

[7] D. Teaney, L. Yan, arXiv:1010.1876| [n ucl-th]]. 

[8] P. Staig, E. Shuryak JaTXiv: 1105.0676] [nucl-th]]. 

[9] E. Shuryak, Phys. Rev. C80, 054908 (2009). see also this proceedings. 

[10] S. A. Voloshin, W. E. Cleland, Phys. Rev. C53, 896-900 (1996). Phys. Rev. C54, 3212-3217 (1996). 

[11] S. Voloshin, R. Lednicky, S. Panitkin and N. Xu, Phys. Rev. Lett. 79, 4766 (1997). 

[12] D. Hardtke, S. A. Voloshin, Phys. Rev. C61, 024905 (2000). 

[13] M. A. Lisa, S. Pratt, R. Soltz, U. Wiedemann, Ann. Rev. Nucl. Part. Sci. 55, 357-402 (2005). 

[14] B. Zhang, C. M. Ko, B. -A. Li, Z. -w. Lin, Phys. Rev. C61, 067901 (2000). 



